Hydration free energies in the FreeSolv database calculated with polarized iterative Hirshfeld charges

Computer simulations of biomolecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in biomolecular systems and are therein described by atomic point charges. In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute’s electron density computed with an implicit solvent model, and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the AM1-BCC and the MBIS atomic charge methods. The latter includes the solvent polarization and presents a root-mean-square error of 2.0 kcal mol–1 for the 613 organic molecules studied. The largest deviation was observed for phosphorus-containing molecules and the molecules with amide, ester and amine functional groups

GTI-space : the space of generalized topological indices

A new extension of the generalized topological indices (GTI) approach is carried out torepresent 'simple' and 'composite' topological indices (TIs) in an unified way. Thisapproach defines a GTI-space from which both simple and composite TIs represent particular subspaces. Accordingly, simple TIs such as Wiener, Balaban, Zagreb, Harary and Randićconnectivity indices are expressed by means of the same GTI representation introduced for composite TIs such as hyper-Wiener, molecular topological index (MTI), Gutman index andreverse MTI. Using GTI-space approach we easily identify mathematical relations between some composite and simple indices, such as the relationship between hyper-Wiener and Wiener index and the relation between MTI and first Zagreb index. The relation of the GTI space with the sub-structural cluster expansion of property/activity is also analysed and some routes for the applications of this approach to QSPR/QSAR are also given

Degeneracy in One-Dimensional Quantum Mechanics: A Case Study

ABSTRACT: In this work we study the isotonic oscillator, V(x) = Ax 2 + Bx −2 , on the whole line −∞ < x < +∞ as an example of a one-dimensional quantum system with energy level degeneracy. A symmetric double-well potential with a finite barrier is introduced to study the behavior of energy pattern between both limit: the harmonic oscillator (i.e., a system without degeneracy) and the isotonic oscillator (i.e., a system with degeneracy)

A semiempirical approach to the intra-phycocyanin and inter-phycocyanin fluorescence resonance energy-transfer pathways in phycobilisomes.

A semiempirical methodology to model the intra-phycocyanin and inter-phycocyanin fluorescence resonance energy-transfer (FRET) pathways in the rods of the phycobilisomes (PBSs) from Fremyella diplosiphon is presented. Using the Forster formulation of FRET and combining experimental data and PM3 calculation of the dipole moments of the aromatic portions of the chromophores, transfer constants between pairs of chromophores in the phycocyanin (PC) structure were obtained. Protein docking of two PC hexamers was used to predict the optimal distance and axial rotation angle for the staked PCs in the PBSs' rods. Using the distance obtained by the docking process, transfer constants between pairs of chromophores belonging to different PC hexamers were calculated as a function of the angle of rotation. We show that six preferential FRET pathways within the PC hexameric ring and 15 pathways between hexamers exist, with transfer constants consistent with experimental results. Protein docking predicted the quaternary structure for PCs in rods with inter-phycocyanin distance of 55.6 A and rotation angle of 20.5 degrees . The inter-phycocyanin FRET constant between chromophores at positions beta(155) is maximized at the rotation angle predicted by docking revealing the crucial role of this specific inter-phycocyanin channel in defining the complete set of FRET pathways in the system

Hydration Free Energies of Organic Molecules in the FreeSolv Database Calculated with Polarized Atom In Molecules Atomic Charges and the GAFF Force Field.

Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges.In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute\u27s electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the MBIS atomic charge method, including the solvent polarization, with a root mean square error of 2.0 kcal mol-1 for the 613 organic molecules studied. The largest deviation was observed for phosphor-containing molecules and the molecules with amide, ester and amine functional groups.</div

Hydration Free Energies of Organic Molecules in the FreeSolv Database Calculated with Polarized Atom In Molecules Atomic Charges and the GAFF Force Field.

<div>Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges.</div><div>In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute's electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the MBIS atomic charge method, including the solvent polarization, with a root mean square error of 2.0 kcal mol^-1 for the 613 organic molecules studied. The largest deviation was observed for phosphor-containing molecules and the molecules with amide, ester and amine functional groups.</div